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Abstract. We have observed the central region of the 
nearby starburst galaxy NGC 4945 with the mid-infrared 
spectrometers SWS and PHT-S aboard ISO. We do not 
find any evidence for the existence of the powerful AGN, 
inferred from hard X-ray observations. The upper limits 
on our AGN tracers [Ne v]14.32/xm&24.3/Ltm and [Nevi] 
7.65/im imply an Ay >160 towards the NLR, assum- 
ing the NLR to be of equal strength as in the Circinus 
galaxy. Other possibilities are discussed. The starburst 
excitation indicators [Nein]15.56/xm/[Nen]12.81/xm and 
Lboi/Liyc suggest that the starburst in the central region 
is at least 5xl0 6 yrs old, and that it accounts for at least 
half of the nuclear bolometric luminosity. The starburst 
might well power the entire bolometric luminosity, but 
the available constraints are also consistent with an up to 
50% contribution of the embedded AGN. With PHT-S, at 
a resolution of «90, we detect strong absorption features 
of water ice, and, for the first time in an external galaxy, 
of CO2 and CO. The same PHT-S spectrum also reveals 
strong emission from the family of PAH features. Finally, 
we have observed and detected several pure rotational and 
ro-vibrational H2 lines, two of which, the (0-0) S(0) & 
S(l) lines, allow us to determine the excitation temper- 
ature (160K) and warm H 2 mass (2.4xl0 7 M Q ). The low 
excitation temperature shows Orion-like shocks not to be 
representative for the entire emission of the central region 
of the galaxy and fairly normal PDRs to be perhaps more 
typical. 
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1. Introduction 

NGC 4945 is a nearby, large (20' x4') spiral galaxy seen 
nearly edge on (i ~ 78°; Ott 1995). At a recession velocity 
of 560km/s it is at the mean radial velocity of the Cen- 
taurus group (Hesser et al. 1984), of which it is believed 



to be a member. Distance es timat es vary between 3.5 and 
4.0M pc (sec Bergman et al. 1992 and Mauersberger et al. 



1996| for discussions). In this paper we w ill adopt a dis- 
tance of 3.9 Mpc (Bergman et al. 1992 ), which implies 
that 1" is equivalent to 18pc. 

NGC 4945 is one of the brightest infrared galax- 
ies in the sky: Sfl2]=24J y, S[ 25]=43Jy, S[60]=588Jy, 



S[100]=1416Jy (Rice et al. [1988Q . The total infrared lu- 
minosity amounts to L(8-1000/xm)=2.95x 10 10 L©, -75% 
of which originates in the central 12" x 9" (Brock et al. 



1988) 



Near infrared observations reveal the nuclear region 
to be the site of a powerful, yet visually obscured, star- 



burst. Br7 (Moorwood et al. 1996a) and Paa (Marconi et 
al. 2000) recombination line maps show the starburst to 
be concentrated in a circumnuclear disk or ring — 200pc 
across (11"). Further evidence for (a period of) strong star 
formation comes from the discovery of a conical structure, 
roughly perpendicular to the galaxy major axis. It is be- 
lieved to be a cavity, vacate d by a starburst- driven s uper- 
wind (Hcckman et al. 1990; Moorwood et al. 



1996a). The 



non-detection of [O ill] within the cone and the absence of 
coronal lines excludes an AGN as the driver of the outflow. 
Clear evidence for the presence of an AGN comes from 



hard X-ray observations (Iwasawa et al. 1993; Guainazzi et 
al. 200C ). The AGN X-ray emission is however heavily ab- 



sorbed by a column density of 10 24 7 cm -2 , which obscures 
the AGN at all optical and infrared wavelengths. Previ- 
ous authors have attributed most of the IR luminosity to 
the starburst (e.g. Moorwood & Oliva 1994; Koornneef 



Israel 1996). Hard X-ray observations with BeppoSAX 



indicate that the bolometric luminosity may as well be 



accounted for by the AGN alone (Guainazzi et al. 2000). 

3cm&6cm ATCA radio ma ps of the central region 
of NGC 4945 (Forbes & Norris |l998|) are dominated by 
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strong nuclear emission, and emission extended along the 
disc of the galaxy. There is also evidence for some fila- 
mentry structure associated with the cavity cleared by the 
starburst superwind. VLBI observations by Sadler et al. 
( |1995| ) reveal the existence of a compact radio core. This, 
as well as the presence of H2O megamasers in a Keplerian 
disc about a ~10 6 M Q black hole (Greenhill et al. |1997| ), 
are taken as further evidence for the presence of an AGN. 

Near infrared observations of molecular hydrogen emis- 
sion in NGC 4945 have been reported by several authors 



Table 1. NGC 4945 results from SWS observations 



1994: Koornneef & Israel 1996; Moorwood et al. 1996a 



over the last 15 years (e.g. Moorwood & Glass 1984; Moor- 
wood & Oliva |1988| Koor nneef Ji993j ; Moorwood & Oliva 



Quillen et al. |1999| ; Marconi et al. |2000p . While fluxes are 
known for eight ro-vibrational transitions accesible from 
the ground (Koornneef & Israel 1996| ) , spatial information 
is available only for the (1-0) S(l) 2.1218/im line. These 
observations show the H 2 emission to be associated with 
the hollow cone, not with the starburst traced in hydrogen 
recombination emission. The absence of a correlation ar- 
gues against photons as the source of excitation. Instead, 
the emission is attributed to shock heating of the molec- 
ular material at the face of the cavity (Moorwood et al. 



19964 Marconi et al. [2000D 



Mid-infrared spectroscopy is much less affected by in- 
tervening extinction than the UV and optical equivalents, 
with A(A)/Ay less than 0.1. Observations of the central 
region of the galaxy, using the mid-infrared spectrometer 
SWS (De Graauw et al. 1996 ) and the spectrophotome- 
ter P HT-S (Lemke et al. |l996| ), both aboard ISO (Kessler 
et al. 1996), are therefore very useful to study the nu- 
clear components otherwise hidden by heavy extinction. 
In Sect. 3.1 we present the results of the search for high 
excitation emission from the AGN. In Sect. 3.2 we study 
the properties of the nuclear starburst. In Sect. 3.3 we 
discuss the dominant nuclear power source. Sect. 3.4 dis- 
cusses the broad emission and absorption features, tracing 
the properties of the interstellar medium in and in front 
of the nucleus. Finally, in Sect. 3.5 we discuss the physical 
conditions and excitation of the warm molecular hydro- 
gen. 

2. Observations 

As part of the Central Program "MPEXGAL" , we have 
observed the central region of NGC 4945 with the Short 
Wavelength Spectrometer (SWS) and the spectropho- 
tometer PHT-S on board ISO. 



2.1. SWS spectroscopy 

SWS grating line profile scans (SWS02 mode) were ob- 
tained on 1996 February 6 for 28 spectral lines in the 
range of 2.42 to 40.34/zm. The spectral resolution in this 
range varies between 11=900 and 2000, corresponding to a 
velocity resolution of 330-150 km/s. Aperture sizes used 
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range between 14"x20" and 20"x33". SWS was ce ntere d 
on the 1.4GHz continuum peak from Abies et al 



( P-987D , 

which coincides with the the position of the L-band peak 
of Moorwood et al. ( |l996a| ). At the time of observation the 
position angle of the major axis of the SWS apertures was 
35.4°, 10° off from that of the galaxy major axis (45°). 

The data reduction was performed using the SWS 
Interactive Analysis system (SIA 



1998 



Lahuis et al 

Wieprecht et al. 1998), putting special emphasis on tools 
to improve cosmic ray spike removal, dark current sub- 
traction and flat fielding. The wavelength calibration of 
SWS is discussed by Valentijn et al. ( 1996a| ). We used cal- 



ibration files as of March 1999. The accuracy of the flux 



calibration is estimated to be 30% (Schaeidt et al. 1996). 
The resulting spectra are shown in Fig [I]. 

In total 17 spectral lines were detected. For another 
11 lines we derived upper limits, using gaussian profiles 
of width equal to other lines of the same (or comparable) 
species, scaled to a peak height corresponding to approx- 
imately 3cr of the noise. Both detections and upper limits 
are presented in Table |l|. 
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Fig. 1. SWS line spectra of NGC 4945. Typical ±lcr error bars are marked. Note that noise increases towards the 
edges of scans 
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2.2. PHT-S spectrophotometry 

We have obtained two low resolution (A/AA ~90) ISO- 
PHT-S spectra, on 1996 October 12 and 1997 August 3, re- 
spectively. ISO-PHT-S comprises two low-resolution grat- 
ing spectrometers covering simultaneously the wavelength 
range 2.47 to 4.87/zm and 5.84 to 11.62/im. The spectrum 
is registered by two linear arrays of 64 Si:Ga detectors 
with a common entrance aperture of 24" x 24". The mea- 
surements were carried out in rectangular chopped mode, 
using a chopper throw of 180". The resulting spectra thus 
are free of contributions from zodiacal light, that would 
otherwise affect the spectrum. The pure on-source inte- 
gration times were 512 and 1024 s. 

The ISO-PHT-S data were reduced using PIAQ 



(Gabriel et al. 1997) version 8.1. Steps in the data reduc- 
tion included: 1) deglitching on ramp level. 2) subdivision 
of ramps in two sections of 32 non destructive read-outs. 
3) ramp fitting to derive signals. 4) masking of bad sig- 
nals by eye-inspection. 5) kappa sigma and min/max clip- 
ping on remaining signal distribution. 6) determination 
of average signal per chopper plateau. 7) masking or cor- 
rection of bad plateaux by eye- inspect ion. 8) background 
subtraction using all but the first four plateaux. 9) finally, 
flux calibration, using the signal dependent spectral re- 
sponse function of Acosta-Pulido ( 1999| ). This spectral re- 
sponse function avoids some deficiencies of the previous 
PIA response function, in particular in the 3/zm region 
near the "Ice" feature. The absolute calibration is accu- 
rate to within 20%. 

The two resulting spectra were obtained at slightly dif- 
ferent position angles about the nucleus. For the first, the 
square aperture was aligned with the galaxy major axis 
(45°). For the second, the position angle was 31.1°. Fig. | 
shows the averaged ISO-PHT-S spectrum. The on-source 
integration times were used as weight factors in the com- 
putation of the average spectrum. 

A number of emission lines can be identified in the ISO- 
PHT-S spectrum. These include 9.66/im H 2 (0-0) S(3), the 
unresolved blend of 6.99/im [Aril] and 6.91/xm H2 (0-0) 
S(5), and 4.05/rai H Bra. For the 9.66/im H 2 (0-0) S(3) 
line, not included in the SWS02 line scans, we measure a 
flux of 5.4xl0~ 20 W/cm 2 , with an uncertainty of 30%. 



3. Results 

3.1. AGN not seen at mid- infrared wavelengths 

NGC 4945 is a peculiar and interesting target for study- 
ing the relation of AGN and star formation in galaxies. 
Clear evidence for hidden AGN activity comes from hard 
X-ray observations. NGC 4945 is amongst the brighest 
hard X-ray emitting galaxies and exhibits variability of 
its 13-200keV flux on timescales of ~10hrs, which clearly 



establishes its AGN origin (Iwasawa et al. 1993; Guainazzi 
et al. |2000| ; Marconi et al. |2000[ ). The AGN X-ray emission 
is heavily absorbed by a column density of 10 24 7 cm -2 
(corresponding to Ay~2600), a high value, but within 
the range observed for Seyfert 2 galaxies (e.g. Risaliti et 
al. 1999). In unified schemes, the X-ray obscuration mea- 
sures a line of sight towards the very center. Obscuration 
towards the NLR probes a different line of sight and is usu- 
ally significantly lower, making the NLR visible in Seyfert 
2 galaxies. 

Mid-infrared high excitation lines are able to penetrate 
a far larger dust obscuration than their optical and UV 
counterparts. They are therefore ideally suited as tracers 
of embedded AGN activity. Mid-infrared emission lines 
like [Nev] 14.32/mr & 24.32^m, [Nevi] 7.65/xm and [Oiv] 
25.9/im are prominently present in the spectrum of all 



Seyferts observed with ISO (Moorwood et al. 1996b; Gen- 
zel et al. 1998). On the other hand, the same emission lines 



1 PIA is a joint development by the ESA Astrophysics Divi- 
sion and the ISO-PHT Consortium 



are also weakly visible towards, for instance, supernova 
remnant RCW 103 (Oliva et al. |1999| ). [Oiv] emission (at 
the few percent level compared to [Ne 11] ) has also been de- 
tected in a sample of starburst galaxies (Lutz et al. 1998), 
again at a level much weaker than seen in typical AGNs. 
The origin of the weak level emission in these sources is 
believed to be shocks. A detection of any of the high exci- 
tation lines discussed above does therefore not automati- 
cally imply the detection of an AGN in NGC 4945. 

We do not detect the lines of [Nev] and [Nevi]. No 
trace of [Ne Vi] is seen in the wing of the nearby PAH emis- 
sion feature (Fig. P . From Fig. [j] it might appear that the 
two [Ne v] lines were indeed detected. However, at the level 
at which the features appear, instrumental effects play a 
significant role. In the 14.35/im line scan a strong fringe 
in the relative spectral response function coincides exactly 
with the expected position for the [Nev] line. Depend- 
ing on the size of the emitting area, the feature may be 
entirely attributed to this instrumental effect. We there- 
fore chose to state an upper limit for the 14.32/im [Ne v] 
line. The feature seen in the other [Nev] scan, centered 
at 24.4/im, was registered by only two detectors, although 
12 detectors scanned over the central wavelength. As vis- 
ible in Fig. [l], the feature is redshifted with respect to the 
NGC 4945 systemic velocity. This shift is not observed for 
any other line we observed. We therefore derive an upper 
limit for this [Nev] line too. 

The only detected high ionization line in NGC 4945 
is the 25.9/im [Oiv] line. An AGN contribution to this 
line is possible — to match the limits on higher excitation 
lines, only part of the [O iv] emission would be related to 
an AGN. The detection of possibly shock-related [O IV] in 
many starbursts (Lutz et al. 1998j ) cautions, however, that 
this may be a more plausible origin of [O iv] in NGC 4945. 
The ratio of 0.033 with respect to [Neii]+0.44x [Nem] is 
above averag e, but well within the range observed for the 
Lutz et al. ( 1998| ) starbursts, also considering that the 
high extinction (Sect. 3.2) will increase the observed ratio 




Fig. 2. The average ISO-PHT-S spectrum of NGC 4945. Upper inset: comparison with the ISO-SWS spectrum of the 
line of sight towards the Galactic center. Lower inset: the ISO-PHT-SL spectrum in linear flux scale 



relative to the intrinsic one. A population of Wolf-Rayet 
stars as origin of the [O iv] emission seems unlikely. Lutz 
et al. ( 199S| ) have shown that the [Oiv] emission would 
have to originate in widely dispersed small Hn regions 
and would have to be relatively strong. [O iv] emission 
at this level has not been observed in local star forming 
regions. A conservative analysis will hence not attribute 
the [Oiv] emission in NGC 4945 to the AGN nor to a 
population of Wolf-Rayet stars. 

The limits on high excitation AGN tracers are con- 
sistent with several scenarios, or perhaps more likely a 
combination of them: 



The Narrow Line Region is extremely ob- 
scured even in the mid-IR. We derive an 
Ay >160 (A(7.65/im)=A(14.3/im)=A(24.3/xm)>4.3) 
to the NLR from a comparison of Circinus and 
NGC 4945 [Nev] and [Nevi] line strengths, under the 
assumption that the galaxies' NLRs are similar. The 
choice for Circinus is motivated in Table ||. 
UV photons from the AGN are absorbed close 
to the nucleus along all lines of sight 
The extreme ultraviolet luminosity of the AGN 
is lower than in Circinus. In comparison to the 
Circinus SED, this would imply a large deficiency in 
UV relative to X-ray flux (Table §). 



Table 2. In X-ray properties NGC 4945 is very similar to 
Circinus. However, IR indicators for a narrow line region 
are missing in NGC 4945 



NGC 4945 



Circinus 



Distance [Mpc] 
H col. dens, to AGN [cm" 
Ay to AGN 
L(2-10keV) [L Q ] 
*/L„(100keV)/Lpi R 
L(8-1000^m) [L ] 
[Nev]14.3^m [W/cm 2 ] 
[Nev]24.3Atm [W/cm 2 ] 
[Nevi]7.65/xm [W/cm 2 ] 



3.9 

5xl0 24 

2600 

78xl0 7c 

0.003 d 

2.2xl0 10 

<7.0xl0~ 

<5.5xl0~ 

<7.3xl0~ 



4 a 

4xl0 24b 
2100 

8.8-44xl0 7b 

0.002 d 

1.2xl0 10a 

440xl0" 21c 

244xl0" 21c 

413xl0 _21c 



Siebenmorgenet al. ( 1997 
Matt et al. (|l999|) 
Guainazzi et al. (200C) 
Marconi et al. (]200C|) 
Moorwood et al. ( |l996b| ) 



3.2. Starburst properties 

Near-infrared broad-band and emission-line imaging has 
revealed the nucleus of NGC 4945 to be the site of a 
sizeable starburst, the presence of which is illustrated 
by the conically shaped starburst superwind-blown cav- 
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ity traced at many near-infrared wavelengths (Moorwood 
et al. 1996a; Marconi et al. 2000| ). Hampered by the large 
extinction even in the near-infrared, age estimates for the 
nuclear starburst are sparse and intrinsically uncertain. 
ISO-SWS offers the possibility for the first time to ob- 
serve the mid-infrared line ratio [Nelll] 15.56/im/[Ne n] 
12.81/im. This ratio, which is much less affected by extinc- 
tion than visible and UV lines, is sensitive to the hardness 
of the stellar radiation field and hence is a good indicator 
for the age of the nuclear starburst. We observed the two 
lines in the same ISO-SWS aperture, which was centered 
on the nucleus (see Table |]). 

To estimate the extinction to the NGC 4945 nu- 
clear starburst we use the ratio of the 18.71/im and 
33.48/im [S in] lines. This ratio is commonly used as a 
density diagnostic for the density range 10 2 ' 5 -10 4 5 cm -3 
and is only mildly dependent on the temperature of 
the emitting gas. Assuming a typical starburst gas den- 
sity of 300 cm -3 (Kunze et al. 1996; Rigopoulou et al. 
1996), the intrinsic ratio should be ~0.43 (i.e. the value 



in the low density limit, computed using the collision 
strengths of Tayal & Gupta 199E). The observed ratio is 
far lower: 0.14±0.06. We hence deduce a screen extinc- 
tion of A(18.7/^m)=1.7lQ'5, w hich, u sing the Galactic cen- 



ter extinction law of Draine ( 1989 ; with A(9.7/im)/E(J- 
K)=0.71), amounts to Ay=36j^. This value is to be con- 
sidered an upper limit in case the [S in] emitting area 
is larger than 14"x27", in which case an aperture effect 
causes the [Sin] ratio to be a lower limit. 

Another independent estimate of the extinction is usu- 
ally obtained from hydrogen recombination line strengths, 
assuming 'case-B' conditions. For NGC 4945 we therefore 
observed the 4.05/im Bra and the 7.46/im Pfa line. Both 
were measured in the same aperture. The ratio Pfa/Bra is 
0.25±0.10, whereas 'case-B' recombination theory predicts 
a ratio of 0.32. The extinction at 7.46/im must therefore 
be similar or slightly larger than at 4.05/im. This indicates 
that the grain composition is unusual and probably more 
similar to the composition found in th e line of sig ht to- 
wards the Galactic center (Lutz et al. 1996 Lutz 1999 ) 
than found towards other parts of our Galaxy. An extinc- 
tion towards the NGC 4945 nuclear starburst can therefore 
at present not be derived from lines in the 4-7/im range. 

The extinction derived for the nuclear starburst is 
somewhat larger than the value we derive for the warm 
molecular hydrogen (see Sect. 3.5). This indicates that 
the warm molecular hydrogen and nuclear starburst emis- 
sion are coming from different nuclear components, the 
latter possibly more enshrouded than the former. With 
the unusual grain composition in mind, it is striking how 
well the Galactic Center extinction law fits our molecular 
hydrogen data, resulting in a smooth excitation diagram, 
even for the H2 (0-0) S(3) line in the center of the 9.7/iin 
silicate feature (see Sect. 3.5). We are therefore confident 
that the extinction correction for the starburst, derived 
using the [S ill] ratio, is also useful. 



In order to determine the excitation of the nuclear 
starburst we apply the extinction correction derived from 
the [S ill] ratio to the observed [Nelll]/[Neii] ratio. The 
extinction corrections amount to A(12.8/im)=1.51 and 
A(15.6/Lim)=1.19. The extinction corrected [Ne III] /[Nell] 



ratio is 0.064±o;o32. Thornley et al. ( |2000D list observed 



[Ne ill] /[Neil] ratios for 26 starburst galaxies, all measured 
in the same ISO-SWS configuration. Clearly, NGC 4945 
is among the lowest excitation targets in their sample. 
Note that the ISO-SWS aperture used is large compared 
to the typical size scales in starbursts. The ratios listed 
by Thornley et al. ( 2000 ) should therefore be regarded as 
aperture averaged. 

For starburst galaxies Lboi/Liyc is another measure 
of the excitation of star clusters. Depending on the up- 
per mass cut-off, the star formation decay time scale and 
the age of the clusters, Thornley et al. ( 2000 ) modeled 
Lboi/Liyc to lie between 3 and 200. The measured values 
for starburst galaxies range between 3 and 50. Below we 
will determine Lboi/Li yc for the NGC 4945 nuclear star- 
burst. We assume Lboi=LiR (i-e. no AGN contribution to 
Lir) and estimate Li yc from the dereddened 4.05/im Bra 
flux. For A(4. 05/im )=1.2 (applying the Galactic center 
law of Draine ( |1989| ) for Ay 



=36±if) 



and Li yc =670 LB ra 

j+26 



we find L lyc =8+^xl0 8 L Q and L bo i/Li yc =28±^. Using the 
12.81/im [Neil] lin e and the empirical scaling Li yc =64 
LNeii (Genzel et al. 1998 ) a similar result is obtained. 

Given the variety of possible star forming histories, it is 
hard to constrain the age of the nuclear starburst (assum- 
ing no AGN contribution to Lir). However, both excita- 
tion diagnostics agree on a low excitation which suggests 
an evolved burst with an age in excess of 5xl0 6 years, 
but would also be consistent with a low IMF upper mass 
cut-off. 

Marconi et al. ( E000D show that it is possible to con- 
struct starburst models which are consistent with their 
near- infrared observations of NGC 4945, but differ by the 
total luminosity generated (their Fig. 4). An instantaneous 
burst would not leave space in the energy budget for a siz- 
able contribution from the hidden AGN, whereas a combi- 
nation of instantaneous burst and constant star formation 
would. We would like to point out here that the latter 
model would be inconsistent with the low [Nem]/[Neii] 
ratio observed by us. Only their instantaneous burst is in 
agreement with both the near-infrared and mid-infrared 
observations. 



3.3. What powers the nucleus of NGC 4945? 

The large extinction towards the nuclear starburst and the 
AGN buried within, makes it very difficult to assess the 
contributions of either component to the nuclear bolomet- 
ric luminosity. 

The optical, near-infrared, mid-infrared and far-infrared 
spectra of NGC 4945 are entirely consistent with 
a starburst-like nature: BLR or NLR high-excitation 
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lines are absent; the starburst excitation indicator 
[Ne ill] /[Neil] has a star burst-like value; the ratios of 6fim 
(ISO-PHT-S), S12 or S25 to S60 or S100 fluxes are all 
very low and consistent with emission from cold dust only. 
Furthermore, the ratio Lboi/Li yc =28^5, is perfectly con- 
sistent with the low excitation of the starburst as deduced 
from [Ne ill] /[Neil]. And last, NGC 4945 lies on the radio- 
far-infrared correlation for starburst galaxies (Forbes & 
Norris 1998[ ). Hence, the starburst might well account for 
the the entire observed bolometric luminosity. 



On the other hand, Guainazzi et al. ( 2000 ), who have 
observed the AGN in NGC 4945 in hard X-rays, compute 
the AGN to be able to account for all the bolometric lumi- 
nosity observed, if it has a typical quasar Lx/Lboi ratio. 
Since there is no such thing as a template AGN spectrum, 
the c onvers ion factor applied, Li_iokev/Lboi ~0.05 (Elvis 
et al. 1994), may have an uncertainty which could easily 



allow for the NGC 4945 starburst to dominate the bolo- 
metric luminosity instead. 

The same uncertanties surround the accretion rate of the 
~1.6xlO 6 M0 black hole inferred from H2O maser obser- 



vations (Greenhill et al. 1997). A high but not implausible 
rate of 50% of the Eddington rate (L Ed d~4.1 x 10 10 L ) 
would suffice to power the observed bolometric luminosity. 
Given the wide range of efficiencies inferred for Seyferts, 
this information does not add anything to identify the 
dominant power source. 

In this complex situation with two potentially dom- 
inant power sources, the most significant constraint on 
their relative weight is the total Lboi/Li yc ratio and its 
implications on the Lboi/Li yc of the starburst compo- 
nent. Ljy C is directly constrained by observations, but 
changes for different assumptions on the starburst and 
AGN contributions to the total bolometric luminosity. If 
there is a significant AGN contribution, (Lboi/Li yc )sb will 
be less than the global value of 28. Values as low as ^3 
which are possible for a zero age massive star population 



with Salpeter IMF (e.g. Leitherer & Heckman 1995) are 



strongly inconsistent with th e low excitation observed in 
NGC 4945. Thornley et al. (|2000|) model [Ne ill] / [Ne 11] 
and Lboi/Li y c ratios of starbursts, taking into account 
clusters of different ages within the ISO-SWS aperture. 
An evolving starburst with [Ne ill] /[Nell] =0.064 as in 
NGC 4945 must have a L bo i/Li yc ^15 (their Fig. 8). This 
limit simply reflects the higher Lboi/Li yc of later type O 
stars and persists if the low excitation is due to an up- 
per mass cut-off rather than evolution. With a lower limit 
of ^15 on (Lboi/Li yc ) s b, the starburst contribution to the 
bolometric luminosity must be at least ^50%. 

We hence conclude that the AGN in NGC 4945 plays a 
secondary although most likely not insignificant role in the 
energetics of this nearby galaxy. Extremely small values 
for the AGN contribution to the bolometric luminosity 
would imply an unrealistically high ratio of Lx/Lb i for 
the AGN. The very low inferred black hole mass, the very 
cold mid- infrared to far- infrared colors, and the absence 



of any clear line of sight towards the AGN, support our 
view that starburst activity dominates AGN activity in 
NGC 4945. 



3.4. Emission and absorption features 

The infrared spectrum of the central region of NGC 4945 
obtained with ISO-PHT-S (see Fig. [2]) presents a new view 
of the ISM in starburst galaxies. Even at the low spectral 
resolving power of R«90 the spectrum is dominated by a 
wealth of emission and absorption features. 

Especially prominent is the family of PAH emission 
features at 3.3, 6.2, 7.7, 8.6 and 11.3^m, which ISO 
confirmed to be common-place in most galactic and ex- 
tragalactic ISM s pectra (e.g. Acosta- Pulido et al 

1999| ; Mattila et al. |l999l Clavel et al 



Rigop oulou et al. 

[2000 ). Nevertheless, the weakness of the 8.6 and 11.3/im 
PAH bands in NGC 4945 is unusual. Consistent with 
Ay~36 and with the strength of the absorption features 
discussed below, we explain this weakness by heavy ex- 
tinction, which will suppress these two features placed in 
the wings of the silicate absorption feature. 

Perhaps the most important result, however, is the rich 
absorption spectrum, indicating that we are observing the 
infrared sources in the central region of NGC 4945 through 
a medium containing molecular ices. Interstellar absorp- 
tions of 4.27,um (2343cm- 1 ) solid C0 2 and 4.68-4.67/zm 
'XCN'+CO are detected, the first time in an extragalac- 
tic source to our knowledge. At our resolving power and 
signal-to-noise we cannot determine the contribution of 
4.62/zm (2165cm" 1 ) XCN, 4.67/^m (2140cm" 1 ) CO ice and 
of gaseous CO absorptions to the 4.58-4.67/im absorption 
complex. The strength of the XCN absorption appears to 
be remarkable, suggestive of ice grain processing in an en- 
ergetic environment (Lacy et al. 1984). We defer a more 
detailed analysis of the XCN/CO feature to a future pa- 
per, which will also include the results of follow-up obser- 
vations with ISAAC at the VLT. A strong silicate feature 
is observed around 9.7^m (see also Moorwood & Glass 
|1984D . A deep minimum is also detected around 3.0/im, 
which is suggestive of water ice (or more precise, the O- 
H stretch) absorption. Table || gives column densities for 
some of the absorption features discussed above. The pres- 
ence and strength of these absorption features is consistent 
with the high starburst obscuration d erived from the emis- 
sion lines (but see also Chiar et al. (200C) for variations 
in the strength of features along lines of sight of similar 
Ay). 

At the resolution of ISO-PHT-S the molecular absorp- 
tion features in NGC 4945 show striking similarities with 
the features seen in the ISO-SWS spectrum of the line 
of sight towards the Galactic center (Lutz et al. 1996; see 
Fig. |^) . Observations at our spectral resolution do however 
not permit a detailed comparison. Regarding the 4.26/im 
CO2 feature it is likely that the feature can be attributed 
to solid state CO2, since high spectral resolution ISO-SWS 
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Table 3. Column densities of solid state features towards the nucleus of NGC 4945 and towards the Galactic center 
(SgrA*). In order to derive the column densities we integrated f t v Av over the width of the band and divided the 
result by the bandstrength A. Nh was determined from Nh = 1-9x10 21 Ay, where Ay=30 for SgrA* and Ay=36 for 
NGC 4945 



Species 


Arest 


A 


^"center 




N 




N/N 


ii 




H 


[cm/mol.] 






[mol./cm 2 ] 












NGC 4945 


SgrA* 


NGC 4945 


SgrA* 


NGC 4945 


SgrA* 


H 2 


3.09 


2xl0 _lba 


0.9 


0.50 b 


2.4xl0 18 


1.3xlO i8b 


3.5X10" 5 


2.3xl0 _b 


C0 2 


4.27 


7.4xl0~ 17c 


0.8 d 


0.72 c 


2.0xl0 17 


1.7xl0 17 ° 


2.9xl0" 6 


3.0xl0~ 6 


XCN 


4.60 

















Hagen & Tie lens (|198l|) 
Chiar et al. (|200i 



Gerakines et al. (11995) 

Determined by fi tting a Gaussian profile followed by rebinning to the ISO-PHT-S instrument resolution 
Gerakines et al. (119991) 



observations of other sources indicate that the contribu- 
tion of gaseous CO2 to the observed feature is very small 
(see in particular van Dishoeck et al. 1996). In the 4.6- 
4.8/im region, the spectra of NGC 4945 and the Galactic 
center differ more strongly, and a more careful inspection 
is required to assess the contributions of gaseous and solid 
CO and XCN. ISO-SWS spectroscopy of the Galactic cen- 



ter (Lutz et al. 1996) clearly shows that what we see at 
low resolution as a relatively shallow and broad feature is 
in fact dominated by individual lines of gaseous CO. Con- 
tributions of a potential underlying solid CO/XCN com- 
ponent are possible but difficult to separate until our high 
resolution follow-up observations have been executed. 



3.5. Molecular hydrogen: physical conditions, excitation 
and mass 

Near infrared observations of molecular hydrogen emis- 
sion in NGC 4945 have been reported by several authors 
over the last 15 years. The most complete setof obser- 
vations is published by Koornneef & Israel (1996), who 
observed 8 ro- vibrational transitions with IRSPEC at the 
ESO NTT. With ISO-SWS and ISO-PHT-S we have ex- 
tended the number of observed lines from 8 to 14 by ob- 
serving the pure rotational transitions S(0), S(l), S(2), 
S(3), S(5) and S(7) as well as the (1-0) Q(3) line. The 
latter was also observed with IRSPEC and can therefore 
be used to determine the proper aperture correction fac- 
tor for the IRSPEC data set. An overview of the observed 
lines is presented in Table [|. 

Information on the spatial extent of the H2 emitting 
region is only available for the 2.12//m (1-0) S(l) line 
(Koornneef & Israel 1996; Moorwood et al. 1996a| Quillen 
et al. |1999| ; M arconi et al. p000[ ). Based on Fig. 3a of Moor- 
wood et al. (1996a) we estimate that more than 90% of 



the (1-0) S(l) emission fits within the smallest ISO-SWS 
aperture (14"x20"). It is not unreasonable to expect the 



H2 emitting area to increase with decreasing H2 temper- 
ature. The aperture sizes used to observe the respective 
H2 transitions increase with increasing sensitivity to lower 
temperature H2. Based on this, we will assume in what 
follows that ISO-SWS and ISO-PHT-S have observed all 
available warm H 2 . Further to this, all three instruments 
were centered on the same nuclear position and viewed the 
nuclear region under more or less similar position angles 
(see Sect. 2). We will use the ratio of the 1-0 Q(3) line 
fluxes measured by ISO-SWS and IRSPEC to scale the 
other IRSPEC lines to the ISO-SWS aperture size. This 
ratio is 2.33. 

From the 14 transitions observed it is possible to com- 
pute the upper level populations for 12 levels. We assumed 
the H2 levels to be optically thin. The excitation diagram 
in Fig. H shows a plot of the natural logarithm of the total 
number of H 2 molecules (N u ), divided by the statistical 
weight (g u ), in the upper level of each transition detected, 
versus the energy of that level (E u /k). The plot shows the 
situation after extinction correction (see below). 

The excitation temperature (T cx ) of the gas is the re- 
ciprocal of the slope of the excitation diagram. If the warm 
H2 is in LTE, the excitation temperature directly corre- 
sponds to the kinetic temperature. As is clearly visible 
from Fig. [| (extinction corrected; see below) the excita- 
tion temperature increases monotonically with upper level 
energy, from 160K for the combination of (0-0) S(0) & S(l) 
to 2200K for the ro-vibrational lines. 

In a highly obscured galaxy like NGC 4945, extinction 
corrections to the H2 emission will be important. The ex- 
tinction can be estimated from the H2 data themselves 
taking into account that any known excitation mecha- 
nism should produce a "smooth" excitation diagram for 
the pure rotational lines, and that transitions originating 
in a common upper level should give consistent results. 
More specifically, we use three criteria: 
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Table 4. NGC 4945 molecular hydrogen data. A(A) is the extinction correction in magnitudes; A u i is the Einstein 
coefficient for the transition from level u to level 1. T u is the upper level energy of level u; g u is the statistical weight 
of level u; N U (^,J) is the number of H 2 molecules in upper level u 



Identification 




r obs 


Instrument 


Aperture 
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N U (V Jl/gf, 
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0.44 


2.94x10" 


ii 


510 


5 


1.27xl0 oii 
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10 
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21 
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0-0 S(2) 
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9 
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9 
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9 
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0-0 S(5) 
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8 
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45 


4.74x10 s7 


1-0 Q(l) 
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IRSPEC 


6x6 


1.72 
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7 
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9 


9.69x10 s6 


1-0 Q(2) 
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5.1xl0~ 21 
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6x6 


1.71 
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7 


6471 


5 


8.47xl0 56 


1-0 S(0) 


2.2235 


4.4xl0~ 21 


IRSPEC 


6x6 


1.98 


2.53x10" 


7 


6471 


5 


1.03x10 s7 


1-0 Q(3) 


2.4237 


3.20xl0~ 20 


SWS 


14x20 


1.70 


2.78x10" 


7 


6951 


21 


5.82 x10 s6 


1-0 Q(3) 
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1.36xl0~ 20 


IRSPEC 
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7 
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21 


5.82 x10 s6 


1-0 S(l) 


2.1218 


1.29xl0~ 20 


IRSPEC 


6x6 


2.15 


3.47x10" 


7 


6951 


21 


5.84x10 s6 


0-0 S(7) 
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l.llxlO -19 


SWS 


14x20 


0.40 


2.00x10" 


7 
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7.13x10 s6 


1-0 Q(4) 


2.4375 


2.2xl0~ 21 


IRSPEC 


6x6 


1.68 
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7 


7584 


9 


2.28 x10 s6 


2-1 S(l) 


2.2477 


2.3xl0~ 21 


IRSPEC 


6x6 
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7 
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21 
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6x6 
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7 
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a Before aperture correction 

b Extinction law 'A' & A v =20 (see text) 

c After aperture correction (IRSPEC data only) and extinction correction. Adopted distance D=3.9Mpc 
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Fig. 3. Excitation diagram for molecular hydrogen in NGC 4945. Different symbols are used to distinguish different 
instrumental origin: diamond: ISO-SWS; square: ISO-PHT-S; and cross: IRSPEC. Results are shown for three different 
dereddening schemes, marked by different shades of grey. Light-grey denotes the combination extinction law 'B' & 
Ay=50 , middle-grey denotes extinction law 'B' & Ay=20 and black extinction law 'A' & Ay=20 
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— The excitation temperature should increase monotoni- 
cally from the lowest to the highest energy levels. This 
sets limits on the extinction correction for the (0-0) 
S(3) line in the center of the 9.7/mi silicate absorption 
feature. 

— The ratio of the (1-0) Q(3) & (1-0) S(l) lines at 2.42 
& 2.12/xm should be its intrinsic ratio determined by 
molecular constants only. The same applies to the (1-0) 
Q(2) & (1-0) S(0) lines at 2.41 & 2.22^m, that origi- 
nate from identical upper levels too. 

— In LTE, the upper level populations normalized by the 
statistical weights should be similar for the 0-0 S(7) & 
1-0 Q(3) lines at 5.51 & 2.42/im, which differ by only 
4% in upper level energy. 

We have varied the extinction and tried several extinc- 
tion laws. We present the most applicable extinction laws 
here: 

— Law A: A(A)cx A -175 for A <8/xm. For A >8 /xm wc 
took the Galactic center law of Draine ( 1989|), with 
A(9.7/mi)/E(J-K)=0.71 (Roche & Aitkcri |198ED and 
E(J-K)=5. 

— Law B: The same as law 'A', except for the range 
A=[2.6,8.8]/xm, where we took the extinction law as 
found towards the Galactic center (Lutz 1999 ). In the 
3-8/im range this reddening law constitutes a signifi- 
cantly higher extinction than usually assumed. 

From Fig. |^ and the criteria defined above, moderate 
extinctions of Ay = 17-23 are clearly preferred. Extinction 
law A provides a somewhat better fit than extinction law 
B. None of the 3 solutions gives a good fit to the (1-0) 
Q(4) data point. In the following analysis, we use the pre- 
ferrred extinction correction of Ay=20^g and extinction 
law A. We note that the extinction to the H2 emitting re- 
gion is slightly less than that to the starburst [H 11] regions 
(Sect. 3.2). This plausibly match es the morphological re- 
sults of Moorwood et al. ( 1996a ), who find the starburst 
in an obscured disk, but the H2 emission extending into a 
less obscured wind blown cavity. 

A rough estimate of the amount of warm molecular hy- 
drogen in the nucleus of NGC 4945 can be derived from the 
level populations of the pure rotational S(0) and S(l) tran- 
sitions. The excitation temperature for the upper levels of 
these transitions (J=2 and J=3) is 160K. Assuming the 
same excitation conditions for the J=0 and J=l levels, we 
compute a warm molecular hydrogen mass of 2.4x 10 7 Mq. 
This is 9% of the total H2 gas mass estimated from CO 
and 0.7% of the dynamical mass interior to the molecular 



Table 5. Warm molecular hydrogen mass estimates for 
the nucleus of NGC 4945 using the best fit power law 
dM/dT=4.43xl0 15 T~ 4 - 793 M /K. The total H 2 gas mass 
estimated from CO amounts to 2.7x10 s Mq 



ring (Bergman et al. 1992; see below). 

As already noted, the excitation temperature changes 
significantly with level energy. This is the consequence of 
the natural fact that the emitting gas will consist of a mix- 
ture of temperatures. The rich NGC 4945 dataset allows 
us to address this in a more quantitative way. Experiments 
with fits assuming a number of discrete temperature com- 
ponents lead us to suggest that a power law might give 



Temp, range 


M(warm H2) 


% of total M(H 2 ) 


[K] 


[M ] 




200-10000 


2.19xl0 b 


0.8% 


150-10000 


6.51xl0 6 


2.4% 


120-10000 


1.52xl0 7 


5.6% 


100-10000 


3.03xl0 7 


11% 


90-10000 


4.52xl0 7 


17% 


80-10000 


7.06xl0 7 


26% 


70-10000 


1.17xl0 8 


43% 


60-10000 


2.10x10 s 


78% 


50-10000 


4.20x10 s 


156% 



a good representation of the mass distribution as a func- 
tion of temperature. We obtain a good fit for the following 
power law: dM/dT=4.43x 10 15 T~ 4 - 793 Mq/K. The qual- 
ity of the fit is shown in Fig. ||. 

Table || gives warm molecular masses for several low 
temperature cut-offs. Since H 2 gas at temperatures below 
70K does not contribute to the (0-0) S(0) flux, nor to 
any of the other pure rotational lines, we cannot verify 
whether our power law mass distribution continues down 
to the lowest temperatures. Nevertheless, we have included 
mass estimates down to a low temperature cut-off of 50K. 
This number is reasonable, since we don't expect the giant 
molecular clouds (GMCs) to be as cold as in the Galactic 
disk (10-20K). Rather we expect conditions as found near 
the Galactic center, where the GMCs are believed to have 



temperatures exceeding 50K (Armstrong & Barrett 1985 ) 
It is interesting to compare our warm molecular hy- 
drogen gas mass estimate with values found in the litera- 



ture (see Moorwood & Oliva 1994 for a review). Bergman 
et al. ( Hgg ), using the inner molecular rotation curve of 
Whiteoak et al. (1990), compute a dynamical mass interior 
to the molecular ring (R<280pc=15.6") of 3.3xlO 9 M . 
The same authors use CO to derive a total H2 gas mass of 
2.7x10 8 Mq for the ring, assuming the rather high kinetic 
gas temperature of 100K. Note that a low temperature 
cut-off of the order 50-60K in our H 2 temperature dis- 
tribution would bring our estimate of the total H 2 gas 
mass in agreement with that derived from the low level 
CO observations. The total H 2 gas mass of 2.7x10 8 Mq 
agrees well with a starburst-like position of NGC 4945 in 
the L m -M(H 2 ) diagram (Young & Scoville [l99l]). 

In Table ^ we list for a number of external galaxies 
and Galactic template sources temperatures and masses 
of the warm molecular hydrogen gas. With a value of 9%, 
NGC 4945 has a warm H 2 gas fraction similar to that 
found for most of the other galaxies listed. Note how- 
ever that the results for NGC 3256, NGC 4038/39 and 
Arp 220 are less well constrained than for NGC 6946 and 
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Fig. 4. Ratio of the fitted and observed H2 excitation diagram for the best fitting power law dM/dT=4.43x 10 15 
T- 4 - 793 M Q /K 



NGC 4945: only for the latter two can the temperature of 
the warm H2 gas be determined from the (0-0) S(0) and 
S(l) transitions directly. For the same reason a compari- 
son of the H 2 gas temperatures is of limited value unless 
they are derived from identical line combinations. Fairly 
low temperatures can be derived from the (0-0) S(0) and 
S(l) lines, 160K and 179K for NGC 4945 and NGC 6946, 
respectively. Limits for other galaxies listed in Table || are 
consistent with a similarly low temperature. This tem- 
perature is well below that observed for an Orion type 
shock (>500K). It is closer to what is observed for the 
same line combination in PDRs (e.g. Orion Bar: 155K, D. 
Rose nthal priv. comm.; S140: 159K, Draine & Bertoldi, 
1999). While a variety of regions may contribute to the 



galaxy-integrated temperature distribution, this compari- 
son clearly shows Orion-like shocks to be not representa- 
tive for the entire emission, and fairly normal PDRs (or 
less energetic shocks) to be perhaps more typical. If ex- 
cited by sho cks (as suggested by the morphology, Moor- 
wood et al. 1996a| ), then the near- infrared H2 emission 
in NGC 4945 must trace a small subcomponent of faster 
shocks. 



4. Conclusions 

The main results of this paper can be summarized as fol- 
lows: 

— The nuclear starburst is heavily obscured by 36l^ 
mag. of visual extinction, as infered from the [S ill] 
18.7/im/33.5/im ratio. 

— The excitation of the nuclear starburst is very 
low, as deduced from excitation indicators 



[Nein]15.56/im/[Neii]12.81/zm and Lboi/Li yc , consis- 
tent with an age of at least 5xl0 6 yrs. Comparison 
with starburst models implies that at least 50% of the 
bolometric luminosity is powered by the starburst. 

— The very low inferred black hole mass, the very cold 
mid-infrared to far-infrared colors, and the absence of 
any free line of sight to the NLR supports the con- 
clusion that the starburst dominates the bolometric 
luminosity. 

— Our mid-infrared ISO spectroscopy does not provide 
any evidence for the existence of an AGN in the nu- 
cleus of NGC 4945. The only high excitation line de- 
tected, the 25.9/im [O iv] line, is most likely produced 
in shocks associated with the nuclear starburst. 

— The AGN, detected in hard X-rays, is unusual in 
not revealing itself at optical, near-infrared and mid- 
infrared wavelengths. Hence, either the NLR is ex- 
tremely obscured (Ay >160), or UV photons from the 
AGN are absorbed close to the nucleus along all lines 
of sight, or the AGN is deficient in UV relative to its 
X-ray flux. 

— Many ISM solid state and molecular features have been 
observed with ISO-PHT-S in the 2.4-11.7^m range. 
Most prominent in emission are the PAH features at 
3.3, 6.2, 7.7 and 11.2/zm. The strongest absorption fea- 
tures are those of water ice, CO2 and CO, seen against 
the nuclear spectrum. These features show striking 
similarities to the absorption features seen towards the 
Galactic center. 

— We have studied the physical conditions, excitation 
and mass of warm H2, combining IRSPEC and ISO 
observations of 14 transitions. We derive a visual ex- 
tinction of 20^3 mag. to the H2 emitting region. From 
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Table 6. Warm molecular hydrogen in external galaxies and Galactic template sources. In column 4 detected lines 
are printed bold, upper limits normal. Toi refers to the excitation temperature computed from the (0-0) S(0) & S(l) 
fluxes. Likewise, T 12 refers to the excitation temperature computed from the (0-0) S(l) & S(2) fluxes. T rot is the best 
fit excitation temperature to several of the lowest rotational transitions. The warm H2 gas mass (M warm ) is computed 
using the gas temperature listed in any of the preceeding three columns. The last column gives the fraction of H2 gas 
in the warm component 



Object 


Type 


Ref 


Rot. lines observed 


Toi 


T12 


T ro t 




% of total H 2 










[K] 


[K] 


[K] 


[M Q ] 




NGC3256 


galaxy 


a 


SO, SI, S2, S5 


150 a 






10 9 


3 


NGC4038/39 


galaxy 


b 


SI, S2 


200 a 






5.6 xlO 7 


8 












405 




8xl0 6 


1 


NGC4945 


galaxy 




SO— S3,S5,S7 


160 






2.4xl0 7 


9 












380 




1.2xl0 6 


0.4 


NGC6946 


galaxy 


c 


S0,S1,S2 


179 b 






5xl0 6 


5-10 


Arp220 


galaxy 


d 


SO, SI, S2, S5 


150 a 






3.5xl0 9 


10 


Orion Bar 


PDR 


e 


S0-S5,S7,S12,S13,S15,S19,S21 


155 










S140 


PDR 


f 


SO— S5,S7,S9,S13 


159 










Orion Peak 1 


shock 


e 


SO, S1-S21, S25 




578 








Cepheus A West 


shock 


g 


SI— S5,S7,S9 






700 






Cepheus A East 


shock 


h 


SO, S1-S8,S9,S10,S11 






740 







a assumed; limits are measured for NGC 3256 (> 140K) and Arp 220 (> 114K) 
b Toi recomputed from the original (0-0) S(0) & S(l) fluxes 



References: a) Rigopoulou et al. (1996); b) Kunze et al. ( 1996 ) ; c) Valentijn et al. (1996b); d) Sturm et al. ( 1996] ); 
Rosenthal (priv. comm.); f) Draine&Bertoldi (|l999[); g) Wright et al. (|l996j); h) van den Ancker et al. (|2OO0|) 
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the (0-0) S(0)& S(l) lines, we compute a warm (160K) 
H2 gas mass of 2.4xlO 7 M , 9% of the total gas mass 
inferred from CO. The excitation diagram is best fitted 
by a power law of the form dM/dT=4.43x 10 15 T -4 793 
Mq/K. The low excitation temperature of 160K shows 
Orion-like shocks not to be representative for the en- 
tire emission, and fairly normal PDRs to be perhaps 
more typical. 
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